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ABSTRACT 

There have been claims in the literature that accretion disks in which a centrifugally driven wind is the 
dominant mode of angular momentum transport are inherently unstable. This issue is considered here by 
applying an equilibrium-curve an alysis to the wind-driving, ambipolar diffusion-dominated, magnetic disk 
model of Wardle & KonigJ ( 119931) . The equilibrium solution curves for this class of models typically exhibit 
two distinct branches. It is argued that only one of these branches represents unstable equilibria and that a real 
disk/wind system likely corresponds to a stable solution. 

Subject headings: accretion, accretion disks - galaxies: jets - ISM: jets and outflows - methods: analytical - 
MHD 



1. INTRODUCTION 

The ubiquity of energetic and highly collimated jets in com- 
pact astronomical objects (ranging from young stellar ob- 
jects to active galactic nuclei) is often interpreted in terms 
of a universal mechanism: hydromagnetic outflows from ac- 
cretion disks (e.g., lLivioll2000l) . In a pioneering paper that 
presented semianaly tic self-similar solution s for centrifugally 
driven outflows, Blandford & Pavne ( 1982) demonstrated that 
such jets could efficiently transport angular momentum from 
the underlying disks and suggested this as an alternative 
mechanism to the radial viscous angular-momentum transport 
that has traditionally been invoked in accre tion-disk models. 
KonigJ lfT989l) and lWardle & Konigll Jl993l hereafter WK93) 
subsequently incorporated the disk into the self-similar wind 
model and suggested that the apparent connection between 
disks and jets may at least in part be a reflection of the fact that 
the vertical magnetic angular-momentum transport is a nec- 
essary ingredient in the accretion process. Additional semi- 
analytic models have since been constructed by several au- 
thors (e.g..lFerreira & Pelletierl 1 993aJbl fl995l OT9951 [79961 
Ferrei rafl997tlOgilvie & Livioll 9981.1200 It ICasse & Ferreiral 
2000a b; Ferreira & Casse 2004) and have served to refine our 
understanding of the equilibrium configurations. However, 
the stability of the derived disk/wind solutions is still being 
debated. 

WK93 and Koni gl & Wardld (119961) argued that wind- 
driving disks should be immune to the most powerful cata- 
loged disk instabilities. In particular, they pointed out that 
disks in which a centrifugally driven outflow transports all 
the liberated angular momentum naturally lie in a stabil- 
ity "window" in which the magnetic field is strong enough 
not to be a ffected by the magnetorotational instability (e.g., 
Balbu s~& Hawlevll9 98) but is not s o strong as to be su bject to 
a radial interchange instability (e.g.. lSpruit et alJl995l) . How- 
ever, iLubow et alJ d 19941 hereafter LPP94) suggested that 
magnetic wind-driving disks may nevertheless be inherently 
unstable. Based on a simplified model, they derived two re- 
lations between the mass outflow rate per unit area and the 
infall speed at a given disk radius: one of these relations 
yielded an S-shaped curve and the other a monotonically in- 
creasing curve. LPP94 argued that disks that lose angular mo- 
mentum smoothly at all radii through a wind that removes 
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only a small fraction of the inflowing mass correspond to the 
two curves intersecting in the middle portion of the S curve 
an d are therefore unstable. This derivation was questioned 
bv lKonigl & WardTel (119961) . who argued that the equilibrium 
model adopted by LPP94 was overly simplified. In particu- 
lar, they applied the LPP94 prescription to the more compre- 
hensive model considered by WK93 and showed that both of 
the resulting equations for the mass outflow rate were inde- 
pendent of the inflow speed and therefore did not yield the 
equilibrium curves invoked in the LPP9 4 argument. 

This issue was recently revisited by iCao & Spruitl ( 120021 
hereafter CS02), who carried out a linear stability analysis 
on an approximate equilibrium disk model. CS02 derived a 
WKB dispersion relation and used it to infer that magnetic 
wind-driving disks are unstable if the wind torque is strong or 
(when the torque is weak) if the rotation is close to Keplerian, 
but that magnetic diffusion stabilizes the disk if this torque is 
small. As interpreted by LPP94 and CS02, the instability re- 
flects the sensitivity of the outflowing mass flux to changes in 
the inclination of the magnetic field at the disk surface: a per- 
turbation increasing the inflow speed causes the poloidal field 
to be bent closer to the disk surface, giving rise to a higher 
mass loss and thereby a higher angular momentum loss; this, 
in turn, allows more mass to flow in, leading to a further in- 
crease in the inflow speed. 

Although the equilibrium configuration examined by CS02 
contains more of the relevant physics than the model adopted 
by LPP94, it is still unclear to what extent the approxima- 
tions made in that work have impacted the final result. For 
example, the field line shape used in determining the location 
of the sonic point and thence the mass outflow rate was not 
calculated self-consistently. This is reflected in the fact that 
its adopted form does not have an inflection point near the 
disk surface — as it should based on the WK93 results. Fur- 
thermore, the normalized magnetic torque at the disk surface 
(r m o in the notation of CS02) is equal to |V, |/2Vk (where V r 
is the mean radial inflow speed and Vk is the Keplerian rota- 
tion speed), which is typically <C 1 even in disks where the 
magnetic torque dominates the angular momentum removal 
(e.g., WK93). On the other hand, the normalized diffusiv- 
ity parameter (?/ in the notation of CS02) is typically 0(1) in 
steady-state disks where the diffusivity balances the inward 
advection of the field by the accreting matter. 1 Since, accord- 
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ing to the CS02 analysis, low-T^o and high-?/ disks are stabi- 
lized by magnetic diffusivity, it is unclear whether one could 
conclude that real disk/wind systems are unstable even if one 
were to accept the CS02 results. 

Numerical simulations offer a potentially effective means 
of resolving the stability question. However, global numer- 
ical studies of a magnetic wind-driving disk incorporating 
a magnetic diffusion mechanism that enables the system to 
reach a steady state have only recently become feasible (e.g., 
ICasse & Re opens 2002, 2004). The results reported to date 
are indicative of stability in that they show that the disk/wind 
system evolves to a near-stationary state, but they clearly do 
not yet settle the issue. In an attempt to shed more light on 
the problem, this paper adopts the equilibrium-curve approach 
first used by LPP94 but applies it to a more comprehensive 
physical model. In particular , using the disk/wind solutions 
of WK93 (§0 see alsoO fl995l fl996l) . it is demonstrated 
(§EJl that not all equilibrium disk/wind structures are unsta- 
ble. It is then argued (§|4j that real disks likely correspond to 
stable configurations. 

2. EQUILIBRIUM DISK/WIND MODEL 

Although the WK93 model involves a number of simplifi- 
cations, it incorporates the full set of structure equations for 
the disk and the wind and thus may be expected to capture 
the basic features of an accretion flow in which the angular 
momentum is removed through a centrifugally driven wind. 
The model assumes a steady-state, axisymmetric and geo- 
metrically thin disk that is in near-Keplerian rotation around 
a central mass. The disk is threaded by a large-scale, open 
magnetic field B that is symmetric about the midplane (where 
B r = B^ = 0, using cylindrical coordinates r, <fi, z). In anticipa- 
tion of an application to protostellar systems, the disk is taken 
to be weakly ionized. In the solutions presented in this paper 
it is assumed that the field is "frozen" into the ion/electron 
component and influences the dominant neutral component 
through ion-neutral collisions (resulting in a relative drift — 
referred to as ambipolar diffusion — between the ions and the 
neutrals). In these solutions it is also assumed that the rela- 
tive drift between ions and electrons (associated with the Hall 
term in Ohm's law) and the radial motion of the field lines 
both vanish. 2 

The disk structure is obtained by solving the mass conserva- 
tion relation as well as the radial, vertical, and azimuthal com- 
ponents of the momentum equation for the neutrals, assuming 
isothermality. The ion number density is determined from 
ionization balance considerations. To further simplify the 
problem, WK93 considered the molecular-gas regime where 
the ion density is constant («j 0.5 cm -3 ); this should typ- 
ically be applicable on protostellar disk scales ~ 10 2 AU, 
where the disk column density is sufficiently low that cosmic 
rays can be assumed to penetrate and ionize the entire vertical 
column. In view of the postulated low ionization the ion mo- 
mentum equation is approximated by equating the ion-neutral 
collisional drag to the Lorentz force. The model also incorpo- 
rates the azimuthal and poloidal components of the induction 
equation as well as the condition V • B = 0. WK93 restricted 

sented by WK93. 

2 WK93 demonstrated that solutions with a nonzero midplane radial ion 
speed Vj ,-o do not differ qualitatively from those in which the field lines do 
not move radially. The actual situation depends on the global structure of the 
star/disk system and may evol ve with time. Note, however, that in the disk 
formation model presented in KrasnoDolskv & Koniel 1 2002) the rotationally 
supported, diffusive disk solution satisfies V; <r o = 0. 



attention to a narrow radial band centered on a given radius r 
and matched the disk solution to a radially self-similar wind 
model, but very similar results are obtained when a global 
self-similar model for both the disk and the wind is consid- 
ered (lLill996l) . 

For specified values of Vk, isothermal sound speed C, 
and midplane mass density po, the disk solution at radius 
r is determined by the values of the neutral-ion coupling 
parameter 77, 3 defined as the ratio of the (Keplerian) dy- 
namical time to the neutral-ion momentum-exchange time 
Tni oc 1 /«;, the magnetic field strength parameter a = Vao /C = 
Bq i '[(4-7T po) 1/2 C] (where B is the field amplitude at z = 0), 
which gives the ratio of the midplane Alfven speed Vao to C, 
and the inflow-speed parameter e = -V r o / C. The solutions de- 
rived by WK93 correspond to the "strong coupling" regime: 
r\ > 1 at all heights. In this case the thermal pressure is not 
much larger than the magnetic pressure at the midplane (i.e., 
a is not <C 1), B r starts to increase already near z = and by 
the time the disk surface (subscript s) is reached it generally 
exceeds \B$. S \ (leading to a strong magnetic squeezing of the 
disk), and the midplan e inflo w speed i s typically of the order 
of C (i.e., e ~ 1). O ^1993 see also lWardld[T997l) showed 
that "weakly coupled" disk configurations, in which 77 < 1 
over the bulk of the vertical column, are also possible. In 
this case a can be <C 1 , B r only starts to grow well above the 
midplane (after 77 comes to exceed 1), with B rs /\B^. S \ usu- 
ally remaining < 1 and magnetic squeezing being relatively 
unimporta nt, and the mass-averaged in flow speed is typically 
<C C. Krasnopolskv & Konigl ( 2002) argued that the latter 
type of disk is likely to arise in the gravitational collapse of 
a rotating, magnetic, molecular cloud core; however, in this 
paper attention remains focused on the strong-coupling solu- 
tions constructed by WK93. 

In the WK93 derivation, the disk solution is matched to 
the "cold" radially s elf- sim ilar centrifugal wind solution of 
Blandford & Pavne ( 1982) by imposing continuity of the 
mass flux and of B, and B^ at the nominal disk surface (the 
height where the disk radial electric field component is equal 
to -VkB z /c, which is required to lie above a density scale 
height in the disk). The details of the transition between the 
diffusive disk interior and the highly conducting wind gas are 
not considered in this model, and neither is the possibility of 
heat injection in a disk corona. 4 Although these issues are rel- 
evant to the determination of the precise range of parameter 
values that give rise to viable solutions (see § |4}, they do not 
affect the basic properties of the equilibrium curves obtained 
in the WK93 model. This is verified, inter alia, by the fact 
that the extension of the WK93 disk solution above the nomi- 
nal dis k surface has similar charac teristics to those of the base 
of thelBlandfor d & Pavnd i ll 9821) wind solution. Any model 
refinements that focus on the transition region between the 
disk and the wind should therefore have no significant effect 
on the stability analysis carried out in §0 

The wind model incorporates the Bernoulli and transfield 
(Grad-Shafranov) equations and its solution yields the flow 
and magnetic field structures above the disk. The wind model 
parameters are k — the normalized mass-to-magnetic flux ra- 
tio, A — the normalized total (particle and magnetic) specific 
angular momentum, and B rs /B z ; however, the condition that 

3 This parameter is the same as fj in CS02. 

4 The implicatio ns of gas heating abo ve the disk surface have been dis- 
cussed by Casse & Ferreiral feOOOtl) and IFerreira & Cassel |2004|) : see also 
IOlmIvk- ,V l.ivi.HI'W4i:n(Ml i 
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FIG. 1. — Mapping of the wind-driving disk solutions (for several values 
of the neutral-ion coupling parameter rf) onto the self-similar wind solution 
space (defined by the values of the mass-loading parameter K and angular- 
momentum parameter A). The field-strength parameter a increases on mov- 
ing counterclockwise along a given curve. 



the wind pass through the Alfven critical point reduces the 
number of independent parameters to two. Since re and A 
can be expressed in terms of the disk solution variables, the 
Alfven constraint serves to reduce the number of parameters 
in the underlying disk model as well. The disk outflow is also 
required to pass through a "sonic" critical point — this yields 
the ratio of the "sonic"-point density to the midplane density, 
which enters into the relationship between the wind and disk 
parameters. 5 

Figure (which reproduces Fig. 11 in WK93) shows the 
mapping between the disk and wind solutions obtained in the 
WK93 model. 6 For given values of rj and C, the disk solutions 
lie along a double-branched curve in the re— A wind parameter 
space. The solutions along each curve are parameterized by 
a, which increases from a value <C 1 at the lower-right end of 
the upper branch to a value ~ 1 at the lower-right end of the 
lower branch. The lower branches of the solution curves end 
on the long dashed curve, below which the outflows remain 
sub-Alfvenic, whereas the upper branches end on the short- 
dashed curve, to the right of which the surface layers of the 
disk are super-Keplerian. As noted in WK93, \B^^\/B Z , and 
hence the magnitude of the angular momentum that the out- 
flow must carry away, increases with decreasing a. Initially, 
as a decreases from ~ 1, B rs /B z (which scales approximately 
as 1/a) increases rapidly, and the corresponding increase in 
the cylindrical radius of the Alfven point (the effective lever 
arm for the back torque exerted by the outflow on the disk, 
which scales as A 1 ' 2 ) increases the value of A and leads to a 
reduction in the ratio of the mass outflow to the mass inflow 
rates (estimated as M oM /M m rj 1/[4(A- 1)]). However, as a 

5 WK93 approximated the outflow at this point as still being largely diffu- 
sive, in which case the relevant "sonic" speed is the thermal speed of sound. 
If, however, the ions and neutrals are sufficiently well coupled to effectively 
form a single fluid at this location, then the critical point corresponds in- 
stead to the slow-magnetosonic point (if th e shape of the field lines at the 
"sonic" point is assumed to be given; e.g., Oeilvie 1997) or the modified 
slow-magnetosonic point (i f the transfield equat ion is also included in the 
calculation — e.g.. lL1199.4lVlahakis eTaTl2 00(». 

6 See Fig. 3 in[L| t_1995f) for the analogous result when the global self- 
similarity solution is extended to also encompass the disk, with the field 
diffusion handled using the generic Ohmic form and a turbulent-diffusivity 
parameterization. 
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FIG. 2. — Disk/wind solutions in the X-a plane. 



continues to decrease, the rate of increase of B r . s /B z declines 
while that of \Bs , S \/B Z increases, and eventually the mass out- 
flow rate must start to increase (with A going down) to keep 
up with the angular momentum removal requirements. The 
transition between these two modes of enhanced angular mo- 
mentum transport: predominantly by the lengthening of the 
lever arm (on the lower branch) vs. mainly by a higher mass- 
loss rate (on the upper branch) occurs at the turning point of 
the solution curve. 

3. STABILITY CONSIDERATIONS 

A turning point in the equilibrium curve typically signals a 
change in the stability properties of the corresponding solu- 
tion. The usefulness of this approach to the disk/wind solu- 
tions shown in Figure ^ is perhaps best demonstrated by re- 
plotting the curves in the A- a parameter space (Fig. 13 based 
on Fig. 12 in WK93). In analogy with thermal stability anal- 
yses, one may think of the parameter a as representing tem- 
perature, the equilibrium value of which is determined by the 
balance between heating and cooling. 7 Within the framework 
of the WK93 formulation, in which the midplane density is 
time independent but the field lines can potentially undergo a 
secular radial drift, the variation of a oc Z?o /^/po is determined 
by that of B z and is governed by the competition between in- 
ward advection and outward diffusion (the analogs of cooling 
and heating, respectively). Specifically, 
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where V; ,• = V r + Vb, r - For the solutions shown in Figures ^ 
and|2 the advection term V r exactly balances the ion-neutral 

7 It is instructive in this regard to compare the present discussion of t he sta- 
bility properties of the WK93 magnetic disk model with Bell & Lin's 1 1994) 
analysis of the thermal ionization instability in viscous protostellar accretion 
disks using the disks' equilibrium vertical structure curves. 
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drift velocity Vd., at z = 0, so the midplane radial ion velocity 
vanishes along the equilibrium curves. 8 

In the case of a self-similar disk that matches onto the 
iBlandford & Pavnel (119821) wind solution, rV; ,B Z scales as 
r" 3 / 4 (e.g., |E3D996) and hence the right-hand side of equa- 
tion Q, when evaluated at z = 0, has the same polarity (posi- 
tive or negative) as Vi. r Q. This implies that the value of a will 
decrease with time if advection dominates diffusion, which 
is consistent with the physical picture of inward gas motion 
bending the field lines and giving rise to a higher value of 
B ns /B z (xl/a. The polarity of V; r o in turn depends on the rel- 
ative magnitudes of the advection (< 0) and diffusion (> 0) 
terms. One can estimate the values of V r and V D r from the 
angular momentum and radial ion momentum equations, re- 
spectively, which can be taken to have a steady-state form if 
the timescale of interest is longer than the dynamical (rota- 
tion) time and if the field is well coupled to the matter (r\ > 1). 
For a thin Keplerian disk with a symmetric field configuration 
one obtains, to a good approximation, 



and 




where h is the disk scale height. As shown by WK93, B r s gen- 
erally exceeds \B^, S \ for a strongly coupled disk, so the con- 
dition of near-hydrostatic vertical equilibrium gives B rs /B z ss 
\[2ja. Further more, \Bj, .JB 7 \ = k(\- 1) from the definition 
of the Blandford & Pavne ( 1982) self-similar wind parame- 
ters. Hence equations 101 and Q imply 

JM-^^A-l). (4) 

VD,rQ 

Equation @ indicates that, if one moves to the right of any 
given equilibrium curve in Figure [2] by increasing A (while 
keeping the other parameter values fixed), then the value of 
IVrol/Vo.rO will go up. In this way we infer that advection 
dominates diffusion to the right of the equilibrium curve (with 
the converse holding to the left of the curve). It follows that, 
if a is perturbed from its equilibrium value, then, by equation 
0, the perturbation will continue to grow if it starts on the 
lower branch of the curve but will be counteracted if it origi- 
nates on the upper branch. From this we deduce that the lower 
branch of an r\ = const equilibrium curve in Figure [2] (or the 
upper branch of the corresponding curve in Fig. Q represents 
unstable solutions, whereas the upper branch in Figure[2](or 
the lower branch in Fig. is stable. 

In the heuristic instability argument given by LPP94 and 
CS02, a decrease in the field inclination to the disk surface 
necessarily leads to an increase in the mass outflow rate from 
the disk. However, as noted in § |2j the equilibrium solutions 
of WK93 have revealed that a decrease in the field inclina- 
tion could also lead to a larger Alfven radius and that these 
two modes of angular momentum transport act in competition 
with each other. As it turns out, the solution branch iden- 
tified here as being unstable is in fact the one along which 
an increase in the angular momentum transport is accom- 
plished through a higher mass outflow rate rather than through 

8 The incorporation of a possible imbalance between field advection and 
diffusion into the WK93 model is analogous to the al lowance for a departure 
from vertical thermal balance in the Bell & Lin 1 1994) viscous disk solutions. 



a lengthening of the effective lever arm. Nevertheless, the in- 
ferred existence of a stable branch indicates that a perturbation 
that reduces the field inclination to the disk does not necessar- 
ily trigger an instability. This may be understood from the fact 
that an increase in B ns /B z also results in greater field-line ten- 
sion, which tends to oppose the inward poloidal field bending. 
Whether a given solution branch is stable or not is determined 
by the extent to which this (as well as any other relevant) 
stabilizing effect can overcome the destabilizing influence of 
increased angular momentum removal brought about by the 
field-line bending. 

4. CONCLUSION 

The stability properties of centrifugal wind-driving accre- 
tion disks are inferred in this paper from a consideration of 
the underlying equilibrium model. Specifically, the ambipo- 
lar diffusion-dominated disk model of WK93 is shown to pos- 
sess equilibrium curves with two distinct branches, and it is 
argued that the turning point between them marks a transition 
between stable and unstable disk/wind configurations. Al- 
though this remains to be verified, a growing perturbation on 
the unstable branch may well correspond to the unstable mode 
identified by CS02 and originally given a physical description 
by LPP94. However, the existence of a generic stable branch 
implies that such disks are not inherently unstable. 

The nonlinear evolution of the instability considered by 
LPP94 and CS02 has not yet been studied, 9 so it is still un- 
clear whether it results in a disruption of the disk. If this were 
the case, then observed systems would naturally correspond to 
the stable branch of the equilibrium solutions. In fact, in some 
cases it may even be possible for a perturbed system on the 
unstable solution branch to undergo a continuous parameter 
evolution into a stable configuration (see Figs. [Jand|2ji. One 
can, however, argue quite independently of the nonlinear out- 
come of the instability that real wind-driving disks are likely 
represented b y stable solutio ns. One argument is based o n 
the results of Ferreira| ill 9971) and lCasse & Ferreiral 112 000a). 
who refined the models used by WK93 and lLil iI1995U199(tI) 
and concluded that the range of viable solutions is reduced; in 
particular, they inferred that all the allowed "cold" solutions 
lie on the equivalents of the lower (i .e., stable) bra nches of 
the curves in Fig ure Q] (see Fig. 3 in Ferreira 1997 and Fig. 
5 in Casse & Ferreira 2000a). Another argument is based on 
the likely formation mechanism of effectively steady disks in 
which the ordered magnetic field originates on large spatial 
scales. In this picture, which naturally applies to the forma- 
tion of protostellar disks from the c ollapse of rotating, mag- 
netic, molecular cloud cores (e.g., Krasnopolskv & Konigl 
2002), the field is advected inward by the accreted matter but 
decouples from the inflowing gas as the ratio of the magnetic 
diffusion time to the inflow time decreases in the vicinity of 
the central mass. Under these circumstances the field lines are 
not strongly bent and B ns /B z is of the order of l. 10 Hence the 
parameter a remains < 1 within the strongly coupled region of 
the accretion flow (B rs /B z would have been ^> 1 had a been 
<C 1) and the corresponding equilibrium configurations lie on 
the upper (i.e., stable) solution branches in Figure|2] 

The equilibrium-curve inferences presented in this paper 
need to be corroborated by an explicit linear stability anal- 
ysis of the WK93 model as well as of more elaborate equi- 

9 It is, however, conceivable that this instabilit y is related to the im plosive 
soliton-like accretion/outflow events discussed by Lovelace et al. 1 1994). 

10 B r ,s/B z = 4/3 in the a symptotic (r — > 0) am bipolar diffusion-dominated 
disk solution of Krasnopolskv & Koniel 1 2002). 
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librium configurations. The arguments in support of sta- 
ble real disks are consistent with the tendency of simulated 
disk/wind systems to approach a near-stationary state (see 
§ 0, although it is worth noting that the equilibrium mod- 
els and disk/wind simulations presented to date have been 
limited to axisymmetric structures. Since the disk and wind 
are potential ly susceptible also to various nonaxisymmetric 
modes (Ye.g.. lKim & Ostrikerl2 000). fully 3D simulations are 
required to conclusively establish their stability properties. 11 
It is also worth keeping in mind that the foregoing conclu- 
sions are based on a consideration of strongly coupled disks, 
whereas (as noted in § |2J the structure of protostellar disks 
may at least in some cases corre spon d to the weak ly cou- 
pled regime discussed bvP lfl996l) and lWardldfi997l) . Since 
in the latter type of disk the bulk of the angular momentum 
is extracted from weakly coupled gas that occupies most of 
the volume whereas the bending of field lines only affects 



strongly coupled gas that resides in thin surface layers, it may 
be expected that such configurations are less susceptible to the 
field-bending instability considered here than their strongly 
coupled counterparts. However, a comprehensive study of the 
stability of weakly coupled disks remains to be performed. 



I am grateful to Z.-Y. Li and R. Rosner for helpful 
comments on the manuscript and to N. Vlahakis for en- 
couragement. This research was supported in part by 
NASA Astrophysics Theory Program grants NAG5-3687 and 
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11 IKrasnopolskv et all <2"p"o"3T) and IQuved etaT] 0553) inferred from 3D 
MHD simulations that steady-state centrifugally driven disk outflows should 
be stable to nonaxisymmetric perturbations. However, these calculations con- 
sidered fixed boundary conditions at the surface of the disk and did not model 
the entire disk/wind system. 
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